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Abstract 

The constraints on the scaling properties of conserved charge densities in the 
vicinity of a zero temperature (T), second-order quantum phase transition 
are studied. We introduce a generalized Wilson ratio, characterizing the non- 
linear response to an external field, H, coupling to any conserved charge, and 
argue that it is a completely universal function of H/T: this is illustrated 
by computations on model systems. We also note implications for transitions 
where the order parameter is a conserved charge (as in a T = ferromagnet- 
paramagnet transition) . 
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I. INTRODUCTION 



There has been some interest in the theory of zero temperature quantum phase tran- 
sitions in condensed matter systems for a few years now 0, particularly in the context of 
metal-insulator transitions @. However, the recent proliferation of experimental systems in 
which such transitions may be observed has lead to a surge in theoretical work. Transitions 
of interest include the superconductor-insulator transition in thin films ||, the transition 
between the plateaus in the quantum hall effect 0], and a variety of magnetic order- disorder 
transitions in the cuprate compounds 0, metal- semiconductor composites ||, and heavy- 
fermion [|7||| systems. 

In this paper we will examine some special properties, in the vicinity of second-order 
quantum transitions, associated with "conserved charges", i.e. observables which com- 
mute with the Hamiltonian. Related issues have been discussed recently by other investiga- 
tors OJTOj , with their focus being on the T = properties of the currents associated with 



conserved charge. We will study here the unusual and remarkable properties of fluctua- 
tions of conserved charges themselves in the finite-temperature quantum- critical |ll|j5| JT2]JT3 
region near the quantum phase transition. 



The quantum-critical region was introduced by Chakravarty et. al 0] in the context 
of the two-dimensional quantum sigma model. An analogous region can in fact be defined 
in the vicinity of any second-order quantum phase transition, as the region where k^T is 
significantly larger than any energy scale which measures deviations of the coupling constants 
from their zero temperature critical values. Note that, somewhat counter- intuitively, the 
quantum-critical region occurs at high temperatures; of course, the temperature cannot 
be so large that it becomes of the order of some high-energy cutoff in the system. At 
short distance/time scales the system displays the scale-invariant properties of some zero 
temperature critical point; at larger scales, the critical fluctuations are quenched by thermal 
effects in a universal manner described in Ref. P,|l3|. Because ksT is large, the thermal 
quenching occurs before the deviations of the couplings from their ground-state critical 
values have had a chance to take effect. Thus, in the quantum-critical region, the dominant 
behavior of the system is described at all scales by the zero temperature critical point and its 
universal response to a finite temperature. Further, the only effect of a finite temperature 
is to impose a finite length h/(kBT) along the imaginary time direction on the quantum 
field-theory of the zero-temperature critical-point; the temperature response of the critical 
point can thus be described by the principles of finite-size scaling [^Tj^|, p!2|,|l3 . 



This paper will examine the non-linear, finite temperature response of a system in the 
quantum-critical region to an external field which couples to a conserved charge. Our moti- 
vation to examine this issue comes primarily from quantum spin systems [ P~2"| , P~5| and heavy- 
fermion alloys although we will attempt to phrase our discussion as generally as pos- 
sible. After some general discussion on conserved charges and their scaling properties in 
Sections |I| and |ITT| , we will present illustrative calculations on a number of model systems 
(Section |TV|). 

Strictly speaking, the considerations of this paper will use only some modest assumptions 
about the zero temperature critical point. In particular we will only require that it be 
gapless with a power-law singularity in the density of low-energy energy excitations. One 
can imagine, particularly in random quantum systems, that this condition may be satisfied 
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even by systems which are strictly not at a scale- invariant critical point. We expect that our 
results will apply to such systems too. However, for definiteness, we will continue to phrase 
our discussion in the language of second-order quantum phase transitions. 

In Section |V|, we will consider a special application of our results to the case where the 
conserved charge is itself the order parameter of the transition: the most familiar example 
of this is the ferromagnet-paramagnet transition in a Fermi liquid. We will show that the 
existing treatment of this transition is fundamentally incomplete, and will indicate the 
restrictions any correct theory must satisfy; we will, however, not provide such a theory 
here. 



II. GENERAL CONSIDERATIONS 

This section will discuss the general constraints that are imposed on correlators of con- 
served charges and currents. These constraints are perhaps most familiar in the particle 



physics context of 'current algebra' ||14|| . We will review these ideas here in a formulation 
designed to address quantum phase transitions in condensed matter systems. Moreover, in 
the latter context, our point of view is different from previous ones flJT0|j, and it therefore 
appears worthwhile to present the complete argument in its full generality. Consider, then, 
the partition function, Z, of the system of interest in the vicinity of the quantum phase 
transition: 

Z = [ V<f) a ex P ( drC[<p a ]) . (2.1) 



<pa(T + L T ) = <f> a {T) V fi 

The Lagrangian £ is a functional of a set of fields <p a which are assumed to be bosonic 
for simplicity - the extension to fermionic fields is straightforward. The fields depend im- 
plicitly on the d spatial co-ordinates x and the imaginary time co-ordinate r. All allowed 
configurations are periodic in r, with period 

L, - if' < 2 ' 2 ' 

where T is the absolute temperature. We will find it more convenient to think of r running 
from — oo to oo with the constraint on the periodicity of the fields, rather than, as is 
conventionally done, restricting attention to the fundamental domain < r < L T . 

Let us now assume that C is invariant (upto a total time derivative), under some 
spacetime- independent symmetry transformation of the fields <p a . In its infinitesimal form, 
this transformation can be written as 

<Pa^<Pa + iVaF^b (2.3) 

where the rj a are the infinitesimal, dimensionless, parameters specifying the transformation 
and the F a are the generators of the Lie algebra associated with the symmetry. These 
generators will satisfy a commutation relation of the form 

[F a ,F?]=tf a ^ (2.4) 
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where the / Q( g 7 are the structure constants of the Lie algebra. 

We now use the usual Noether argument to identify the charges and currents. Make the 
transformation ( |2.3| ) on C, but with the r\ a spacetime- dependent. In general, any variation 
in the action under this transformation, will depend, to linear order, on the derivatives of 
the 7] a . We therefore have for small r] a (and again upto a total time derivative) 

£ — > C + ih J d d xd v j] a q m {x, r) (2.5) 

where the index \i extends over the d + 1 spacetime co-ordinates. The co-efficients, q Ta 
of d T rj a are, of course, the conserved charge densities associated with the symmetry under 
consideration, while the q xa are the associated currents. (At this point, it is conventional in 
some field theory books to identify the q m with (j) a F^ b 5C/5(d^(pb); we caution the reader that 
this latter form fails for the coherent state path integral of quantum spins - the definition 
( |2.5|) is more generally valid.) 

We are interested in the special constraints that apply to correlation functions of the 
q m . To this end, we place the system in external fields A^ a which couple to the q m ] the 
correlation functions can then be obtained by taking appropriate functional derivatives w.r.t 
the A^. While it is sufficient to simply add a linear coupling A^ a q^ a to C to achieve this, 
we will find that this approach is not the most convenient in deriving the Ward identities. 
The following approach is found to be the most direct: Generalize the Lagrangian C[4> a ] to 
the field-dependent C[4> a ,A^ a ] and evaluate 

Z(A tMX )= [ D0 a exp(-~ [dTC[</) a} A^\). (2.6) 

The new C[4> a , A^] is chosen such that it is invariant (upto total time derivates) under 
spacetime- dependent transformations of the form (|2.3|) accompanied by the following trans- 
formation of the A m 

A^a ► A^ a + d^rj a — fafi^pA^ (2-7) 

In other words, we have promoted the global symmetry to a gauge symmetry, and the A^ a 
are the non-abelian gauge connections. 

Let us now examine a few simple examples of the above construction. 



1. Non-relativistic electrons in a magnetic field 

Non-relativistic spin-1/2 electrons, c a (x,r) (a =|, j), i n an external magnetic field H 
rjafrddy/dxp (a, /3,7 = 1,2,3) are described by the following Lagrangian 



C 



d d x 



h jdCa_9Haj a Ch _^_ 

hCa dr 2 Ca ° ahCh 2m 



_d_ 

3Xr 



%Qj a I C a 



+ 



(2.8) 



where the ellipses indicate terms without any derivatives, the a a are Pauli matrices, and 
g is the gyromagnetic coupling. The paramagnetic and diamagnetic couplings of the field 
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to the electrons play totally distinct roles in the present symmetry analysis. There are two 
distinct conserved charges - the total number and total spin of the electrons. The first is 
associated with the U(l) symmetry 

c a ^c a + ir}c a . (2.9) 

Upon gauging this symmetry we see that the a a are the spatial components of the fields 
introduced above 



a a + ^- (2 . 10 ) 



The second is the spin-rotation symmetry 



c a -Ca + iy^* (2.H) 
which when gauged leads to the transformation 

gH a ^> gH a - -d T r] a - gEapyiipHy, (2.12) 

on the magnetic field (e is the totally antisymmetric tensor). Thus igH a /h is the r- 
component of the non-abelian SU{2) gauge field, A m , associated with the SU(2) spin- 
rotation invariance. Note that A Ta is purely imaginary. We shall mainly focus on the 
consequences of this second symmetry in this paper. 



2. 0(3) sigma model 



This model is a popular long-wavelength description of low-lying spin excitations in an 
insulating antiferro magnet. In the presence of an external magnetic field H a , the Lagrangian 



takes the form 15 



2g 



— ( d T n a - —e aab H a n b ) + (d x n a 



(2.13) 



where n a is a 3-component, real, unit-vector representing the local orientation of the anti- 
ferromagnetic order parameter. £ is invariant under an 0(3) symmetry under which 



(2.14) 



while H a continues to transform as in (|2.12|) and thus igH a /h is the r-component of a 0(3) 
non-abelian gauge field. 
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3. Quantum spins 



The symmetry analysis of the path-integral of quantum spin systems is somewhat more 
subtle, but our general discussion has been phrased carefully to include this case. As was 



first shown by Haldane [TfJ, the path integral of any quantum spin Hamiltonian involves the 
Lagrangian 

£ = thSj2W b (Q aj )^- + H(Q aj ) (2.15) 

where S is the half-integral/integral magnitude of the spin, the Q a j are unit 3- vectors on sites 
j representing the instantaneous orientation of the spin, and W a is any function satisfying 

Sabc^rW c = Q a (2.16) 
o\ l b 

(we have momentarily dropped the site index j). The Hamiltonian TC does not involve any 
time derivatives, and is spin-rotation invariant. Let us now make a space- independent, but 
time-dependent rotation of all the spins 

Q a -> tt a - r] c (r)e cab n b , (2.17) 



Inserting this into C, using (|2.16 ) and the unit-length constraint on Q, simple manipulations 



show that, upto a total time derivative, 

C^C + ihS^ -Yn ai (2.18) 
dr ^ 

By our prescription ( |2.5| ) this identifies S J2j ^aj as the conserved total spin. In the presence 
of a magnetic field the action is clearly 

c = ihs^w b {n aj )^- g H a Y,n aj + H(n ad ) (2.19) 

3 i 

This is now invariant under time- dependent gauge transformations with H a transforming as 
in ( 12.121) . Note that the 'rule' of replacing derivatives with covariant derivatives does not 



hold in this case - our formulation is however still valid. 

We now return to the general considerations. We will consider first the Ward identities 
satisfied by correlators of the conserved charges and currents. This will be followed by a 
discussion of properties properties of the system in a time-independent external field. 

A. Ward Identities 

An important property of the functional Z(A^ a ) in (|2.6|) is that 

Z(A m ) = Z{A m + d^ a - f a ^np A m) ( 2 - 20 ) 
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for any spacetime dependent gauge transformation rj a such that (|2.3| ) is consistent with the 
boundary conditions a (r + L T ) = 4> a { T )- (This follows from performing (|2.3| ) on the a 
dummy variables of integration, followed by ( |2.7| ), which leaves the action invariant.) We 
expand (|2.2CI|) to linear order in t] and obtain the key Ward identity 

dA^aix.T) 5A^{X,T) 

The left-hand-side of this equation is simply the divergence of the conserved charge and 
currents. The right-hand side is the analog of the 'streaming' or 'Poisson-bracket' terms [T7| 
in the theory of the dynamics of classical phase transitions; this term dictates that the 
conserved charge undergoes a uniform precession under the presence of the external field. 

In the following we will mostly be interested in constraints on correlators of the conserved 
charges q Ta under conditions in which only the r component of the A m is non-zero. By 
integrating (|2.21|) over all space we can obtain a constraint on these correlators 

id o ^(^to) A xa = 0) r d 5Z(A TCe , A xa = 0) 



9 ' \Z£r) = I d X IMX ' T) SZ** T) (2 ' 22) 

A particularly useful consequence of ( 2.22Q is the constraints it places on two and three-point 
functions of the conserved charges. If we make the expansion (restricting, for simplicity, to 
a translationally invariant system) 

T(A Ta ,A xa = 0) = J d d qdujG(q,u)A Ta ( y q,uj)A Ta ( y -q, -u) + 

J d d q 1 d d q 2 du} 1 du} 2 T al3 ~ / (q 1 , q 2 , UJi, uj 2 )A Ta (qi, LUi)A T p(q 2 , u 2 )A T1 (-q 1 - q 2 , -oj x - u 2 ) 

+ ••• (2.23) 

where the qi and u>i are momenta and frequencies and T is the free energy density, we see 
from (|2T22|) that 

+ u 2 )T a ^(q, -q, u x , u 2 ) = f a ^(G(q, u x ) - G(q, u 2 )) (2.24) 
This identity will be useful to us later in our study of ferromagnets. 



B. Time-independent, uniform, external field 

We will consider explicitly only the case of a time- independent, uniform, A Ta field; the 
spatial components A xa will be taken to be zero. As was clear from the examples considered 
above, the A ra corresponds to an imaginary external magnetic field in spin systems. To 
emphasize this we will use the notation 

A ra = (2.25) 



As in Section pTA] we attempt to 'gauge away' the field dependence of Z(H a ) = Z(A Ta , A 



xa 



0) for the case of a time-independent H a . From fl2.7| ) it appears that we should choose 
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(a generalized Josephson equation). However the corresponding transformation (|2.3| ) on the 
<pa necessarily modifies the boundary conditions. We have therefore 

Z{H a ) = { P0 a exp (~ / drC^a}) . (2.27) 

Thus the sole effect of the field H a is to put a twist in the periodic boundary conditions on 
4> by an imaginary angle igHL T /h. 



III. SCALING PROPERTIES NEAR QUANTUM PHASE TRANSITIONS 

We will focus almost all our subsequent attention in the 'quantum-critical re- 
gion' Hliy5il2lJl3l where /cgT is much greater than any intrinsic low-energy scale associated 



with the deviation of the ground state from criticality (we must of course not make fcgT so 
large that it becomes comparable to ultraviolet cutoff's in the system). In this region, the 
leading T dependence of all observables is specified by properties of the T = critical point. 
In the following, we will therefore neglect the deviation of the ground state from criticality, 
although the extension to including its consequences are quite straightforward. Also, we will 
mostly consider the case of a uniform, time- independent field A TOl ^ 0, A xa = 0, and refer 
to the external field using ( |2.25| ). 



We consider the properties of the the free-energy density T = —(h/(L T V))\ogZ = 
— (ksT/V) log Z (V is the spatial volume of the system which is assumed to be infinite) 
as a function of T and H. Consider first the case H = 0, and T close to 0. The only 
effect of a finite T is in the imposition of a periodicity in <fi with period L T on the critical, 
scale-invariant theory at T = 0, H = 0. The hypothesis of finite-size scaling [18| predicts 
the following temperature dependence in T 

TiT, H = 0) = JF(0, 0) — c\T p (3.1) 

There is no general expression for the exponent p, or the constant c\. However, if the system 



is below its upper critical dimension, the hyperscaling hypothesis |18| states that the scaling 



dimension of T is identical to its naive engineering dimension: this yields 

d 

p=l + - (3.2) 

z 

The 1 contribution is due to the 1/L T prefactor in the definition of J 7 , and the remaining 
d/z contribution is from the 1/V. The dynamic-critical exponent z expresses the anisotropic 
scaling between space and time directions. The pre-factor c% in ( |3.1| ) is in general non- 
universal. For the special case of a relativistic field theory we have z = 1 and the c a becomes 
universally related to the velocity of the low-lying excitations; in this case in d = 1 the 
number C\ is closely related to the central charge of the conformal field theory describing 
the critical point. 
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Now consider the effect of a time-independent external field H. From ( |2.27|) the only 
effect of if is a twist in the r boundary conditions on the system. The excitations responding 
to the change in the boundary conditions will be precisely the same low-energy modes which 
led to a size (T) dependence of T in ( |3.1|) : the only effect of a finite H should therefore 
be a modification of the term proportional to c\ in ( |3.1| ). Furthermore, as the twist is a 
long-wavelength effect, the modification of c\ should be 'universal' (i.e. independent of all 
microscopic details) function of the 'angle' of the twist igHL T /h. Alternatively, this is simply 
the statement that all finite-size scaling corrections are universal functions of 'geometrical' 
properties of the sample like aspect ratios, shape, nature of boundary conditions etc. The 
fact that the angle of the twist is imaginary should not be too disturbing - the process of 
analytic continuation commutes with all scaling arguments, and one can just lift the scaling 
forms from those of real twists. We have therefore 

T(T y H) = F(0, 0) - c 2 T*Q (0;) (3.3) 

where c\ = C2^(0). The value of Q(0) will be chosen at our convenience, but the function Q(r) 
is otherwise universal (we use r = gH/iksT) below). Note in particular that the argument 
of the scaling function is precisely gH/ksT and there are no arbitrary scale factors in the 
argument. There is no guarantee that the function Q(r) is analytic for finite, positive values 
of r. In particular, some systems may undergo a phase transition at a finite H, which will 
then correspond to a (universal) singularity in fi(r); we will see an example of this in the 
model calculations below. 

The form of Q3.3|) implies immediately that the scaling dimension of H (or equivalently 
A TOl is precisely the same as that of T. In other words, under a scaling transformation which 
rescales spatial lengths by a factor s 

A' Ta (x',r') = s z A Ta (x,T) (3.4) 

where x' = x/s and r' = r/s z . Exactly parallel arguments can be made for the spatial 
components of the A m by thinking about the properties of the system in a geometry which 
is finite in the spatial directions, but infinite along the time direction - this will yield the 
scaling dimension of A xa : 

A' xa (x',r') = sA xa (x,T) (3.5) 



We emphasize that that none of the results ( |3.3j ), ( |3.4| ) or (|3.5|) rely upon the validity of 
hyperscaling. 

In the presence of hyperscaling, one can go further, and also deduce the scaling dimen- 
sions of the conserved charges and currents. The q m and the A m are conjugate variables 
and their product should therefore have the same scaling dimension as the free energy (which 
is z + d) . We have therefore 

9™ fa', r ') = s d q Ta {x, r) q' xa (x', r') = s d+z - l q xa (x, r) (3.6) 

only if hyperscaling is valid. 
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A number of strong experimental consequences now follow from ( |3.3|) . We can imme- 
diately obtain scaling forms for the 'magnetization' M = —OT/dH and the specific heat 
C v = -Td 2 F/dT 2 : 

where the universal functions Qm,^c are both simply related to linear combinations of Q 
and its derivatives. Notice that it is the same non-universal number c<i which appears in 
both M/H and Cy, and there are no other non- universal quantities; the only choice that 
had to be made was in the value of Q(0). All dependence on this choice, and hence C2, can 
be eliminated by considering the dimensionless generalized Wilson ratio, W 

Ws &m =thr (m ( 3.8) 



9 2 C v " \k B T 

which is a fully universal function of H/T. We emphasize that the universality of W did 
not rely on hyperscaling. Experimental measurements of this ratio can thus provide us with 
strong tests of various theoretical scenarios, and also determine if different experimental 
systems are in the same universality class. We note that the universality of the Wilson 
ratio as H —>■ has also been noted recently for the incremental thermodynamic response 
of impurities in Fermi liquids [|B5|: these models map onto boundary critical phenomena, 



whereas we have been considering the bulk response of a macroscopic critical system. 



IV. MODEL CALCULATIONS 

We will now illustrate the general principles described above by model calculations on 
a number of systems. We begin with the simplet realization in the theory of Luttinger 
liquids; in this case the function fiw(r) will turn out to be independent of r. None of 
the remaining models will have this property. We follow this by a second simple system 
- a dilute fermi gas - which also satisfies the scaling ansatzes. We will then examine a 
simple phenomenological model of a very complicated system - the Bhatt-Lee [[HJ model 



of random quantum antiferromagnets. Finally we will present a self-contained analysis of 
a truly interacting system: the O(N) sigma model, whose main applicability is to the low- 
energy properties of clean, quantum antiferromagnets. 



A. Luttinger Liquids 

We begin by presenting the simplest illustration of our results in the Luttinger liquid 
theory of the low temperature properties of a dense one- dimensional gas of spin-1/2 fermions. 
In this case we are considering a whole critical phase, rather than a critical point. 

The low energy action of the Luttinger liquid can be expressed in terms of two dimen- 
sionless scalar fields, 9 p , 6 a 



C 



h 

2^ 



dx 



k p u p (d x e p y + —(d T e p y + M Ua (dj a 



a 



u. 



(4.1) 
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where lip, U(j £1X6 the charge and spin excitation velocities, and K p , K a are dimensionless 
couplings which determine the exponents of the Luttinger liquid; we have used here the 
notation of Ref. [^TJ. Spin rotation invariance requires K a = 1. 



In the presence of an external magnetic field couping via the Zeeman term to the spin- 
1/2 fermions, C gets modified by the replacement d T 6 a — > d T 9 a — igH/(^/2h). Computing 
the action of the free field theory £ at finite temperature is now completely straightforward. 
We get 

F(a > T) = -^it {,jH)2 + ¥ s / £ (** H + <«?) + ** H + «?* 2 )) < 4 - 2 > 

Note that the H dependence of T is rather simple and has decoupled completely from its T 
dependence: this is a special feature of the present model. The frequency summations and 
integrals can be performed exactly and yield a result consistent with ( |3.3|) which is: 

^ W D = ^(o,o)-^n(i|) 

The scaling properties of the magnetization and the specific heat now follow. In particular, 
we obtain for the generalized Wilson ratio 

°"M ^ ^HW (4 ' 4) 

As stated above, Qw is in fact independent of r. This is a special feature of the Lut- 
tinger /Fermi liquid that does not generalize. The Wilson ratio has most often been consid- 
ered in the past in the context of Luttinger/Fermi liquids, and this is perhaps the reason why 
its universal, non-trivial, dependence on the ratio H/T at generic quantum-critical points 
has not heretofore been pointed out. 



B. Dilute Fermi Gas 

Consider a gas of fermions (with spin j) in d dimensions described by the following 
Hamiltonian 

where c& annihilates fermions with momentum k, and Ti,^ contains only repulsive interac- 
tions. This model has a T = quantum phase transition as a function of /i at /x = 0. The 
density of fermions vanishes for fi < 0, and increases as ~ fi d ^ 2 for fi > 0. The scaling prop- 
erties of this quantum transition are very similar to those of the corresponding transition for 



bosons which has been studied elsewhere [22]. From this analysis [E2[ we may conclude that 



the exponent z — 2. Also, it can be shown that the interactions in TC[ n i are irrelevant at 
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this transition for d > 2 (they are infact also irrelevant below d = 2 for spinless electrons). 
Thus, for d > 2, we may compute the scaling properties of the free energy in the free fermion 
model. 

The conserved charge we focus on here is the density of the fermions. The field conjugate 
to this density is fi and therefore plays the role here of the 'magnetic' field. Thus consistent 
with ( |3.3|) the free fermion free energy density obeys 

h^t) = -{23 + m B rr^ (frf 7 ^ 

Q{t) = f JLL. log (i + e ~v 2 +r\ ( 4 6 ) 



(2n) d 

Unlike Section [LV A[ note that O(r) is quite a non-trivial function of r, and leads to corre- 



spondingly non-trivial r-dependences in the scaling results for the density (which plays the 
role of 'magnetization'), specific heat, and Wilson ratio. 



C. Bhatt-Lee model. 



This is a simple phenomenological model of the spin-fluid phase (i.e. no spin-glass or- 
der) of spin-1/2 random antiferromagnetic spin systems |20[| . It has been quite successful 
in describing experiments in lightly doped semiconductors [p3 . We now show that this 
model in fact satisfies all of the constraints discussed above on quantum-critical spin fluctu- 
ations. Thus the entire spin-fluid phase may in fact be critical in random systems, and not 
just its transition to a magnetically ordered state. Additional evidence for such a scenario 
has appeared in recent solutions of random Heisenberg antiferromagnets with infinite-range 
interactions . 

The Bhatt-Lee model describes the random antiferromagnet as independent pairs of 
spins which have an antiferromagnetic exchange interaction J with probability P(J) ~ J~ a . 
The exponent a is estimated from numerical work to be approximately 0.6 in d = 3. The 
free energy of this model in an external field H, is obtained by summing the contributions 
of each pair of spins and is therefore 



T = To-k B T / dJP(J) log 



1 + e 



-J/k B T, 



1 + 2 cosh(g H/k B T) 



(4.7) 



This can easily be collapsed into the scaling form (|3.3j ) with p = 2 — a and the universal 
scaling function Q(r): 



n(r) 



dyy~ a log l + e~ 2/ (l + 2coshr) 



(4.8) 



The value of Q(r = 0) has been chosen for convenience; apart from this single scale, the 
function Q(r) is otherwise universal. If we assume hyperscaling then we get the dynamic 
exponent 



d 

1 — a 



(4.9) 
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The integral in (|4.8|) cannot be evaluated exactly, but we quote some useful asymptotic 
limits: 



«m = - r{ l-?^H: 3) r - (4.10) 

w [ r z /((2 — a)(l — a)) + 7rr a /6 r — > oo v y 

where Li p (2;) is the polylogarithm function, defined by analytic continuation of the series: 

00 n 

Li P (-) = £ ~ P (4-11) 

n=l '* 

The scaling functions for the magnetization (fijf), specific heat (Oc), and the Wilson ratio 
can now be easily obtained by taking suitable derivates of Q(r): The results are plotted 
in Figs |l] and ||] for the value a = 0.6 (a function closely related to Qm was evaluated and 
compared with experiments in Ref. f23|). As we noted earlier, the results for Q\y are totally 
independent of any choice of an overall scale - we state below the asymptotic limits of flw- 

2Iii_„(-3) 

r — > 



n w (r) = 3 ( X - a ^ 2 - « )Li2 - Q( " 3) ' (4.12) 

~2?i V r ^ 00 

Note also in Fig |2] the non-monotonic behavior of Q\y between these two limits. 

Our identification of the exponent z also allows us to make a new prediction on the 

temperature dependence of the spin diffusion constant D. The scaling dimension of D is 
z — 2, leading to the low temperature dependence 

D ^ rpl-2/z _ rp(d-2 + 2a)/d^ (4.13) 



D. O(N) sigma model in 2+1 dimensions 

We first generalize the 0(3) sigma model of Section |11 2j to the O(N) model by allowing 
n a to have iV components, a = 1 . . . N. The external field H must now generate one of the 
rotations of the O(N) group. These rotations can be built out of combinations of rotations 
in the N(N — l)/2 different hyperplanes in iV dimensions. For general N, unlike N = 3, not 
all such rotations are equivalent, and cannot be transformed into each other by a change of 
co-ordinates: this is related to the presence of more than a single Casimir invariant in the 
O(N) group. We will therefore choose a specific orientation of the magnetic field to facilitate 
a simple large iV limit: other orientations of the magnetic field will have physically different 
properties (for N > 3). We choose a magnetic field to generate a simultaneous rotation 
by the same angle in the (1, 2), (3, 4), . . . ((2p — 1)N, 2pN) hyperplanes, with no rotation in 
the remaining N(l — 2p) hyperplanes; the fraction p is chosen such that pN is an integer, 
and p < 1/2. The large N limit will be taken with p fixed. Clearly, the model relevant to 
collinear quantum antiferromagnets is iV = 3, p = 1/3. These considerations lead to the 
following action for the O(N) sigma model in a magnetic field: 
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N f f N 2 1 PN 2 2 

S = — dr d 2 x (Vx\) 2 + — 1(^20-1 ~~ igHn 2a f + {d T n 2a + igHn 2a -if 

^ J J a=l C a=l 



1 V 



+3 £ ( 9 - n «) 



C a=2pN+l 



(4.14) 



where the coupling constant t determines the strength of the quantum fluctuations and c is 
the spin-wave velocity. 

The T = phase diagram of S can be deduced by a straightforward extension of the 



methods of Ref. [[Kj and the d = 1 analysis of Affleck |25]; the results are summarized 
in Fig |3|. In zero external field, there are quantum disordered and Neel ordered phases 
separated by a critical point at t — t c . This critical point has z — 1. The quantum 
disordered phase has a gap A which vanishes near t c as A ~ (t — t c ) v . For finite H the 
Neel ordered phase transforms into a second ordered phase in which the spin-condensate 
is preferentially oriented in the pN planes in which the field generates rotation. For the 
physical case pN = 1 this phase has XY order and is so identified in Fig |[ The transition 
between the finite H ordered phase and the quantum disordered phase occurs exactly at the 
field H = H c where the zero-field gap A equals gH. This quantum transition has z = 2 and 



is studied in some detail in a separate paper [26 



The T 0, H properties of S are quite different for the cases pN = 1 and pN > 1. 
We discuss first the case pN = 1, which is summarized in Figs |]-|6|. The T = XY order 
survives at finite temperature as quasi-long-range order. There is a Kosterlitz-Thouless 
transition from this state at a temperature T KT to a fully disordered state. The dependence 
of Tkt on H depends crucially on the value of t. We found (See Figs ^-^) 

2np s /log{p s /H) t<t c 

KbTkt = { g{H _ Hc) iog(A/(lf - H c )) ' ' (415) 



loglog(A/(F-F c )) 



t > t r 



The result for t < t c can be deduced from the results of Nelson and Pelcovits |?7] on a 
closely related model; here p s is the fully renormalized spin stiffness of the ordered state of 
the T = 0, H = 0, sigma model. The situation for t > t c follows from the work of Popov |28 



and is discussed in more detail elsewhere [2^]. Our main focus here is on the t = t c case: the 
finite T, finite H properties can then be deduced by applying the scaling methods of this 
paper to the z = 1 critical point at t — t c . The free-energy of the model continues to satisfy 
( p.3[) . The existence of a finite T Kosterlitz-Thouless transition implies that the function 
fi(r) must be non-analytic at, say, r — K: this leads to the result above for Tkt at t — t c . 
Moreover as there are no non-universal factors in the scale of r, the number K, is universal. 

The finite T properties for pN > 1 are simpler - there is no phase transition at any finite 
T. The windows of quantum-critical behavior with z = 1 (as in Fig |5|) and z = 2 (as in 
Fig are however still defined. 

We will now present the N = oo computation of the universal function Q(r) in the 
vicinity of the T = 0, H = critical point at t = t c . As the large limit is taken with 
p fixed, we necessarily have pN > 1 and there is no finite temperature Kosterlitz Thouless 



14 



transition as in Fig |5|. Nevertheless, from the insight gained in Ref . |H| , we expect that our 
results for Q(r) are reasonable accurate for the physical case iV = 3, p = 1 provided r 3> 1 
or gH 3> ksT. 

The technical steps in obtaining the N = oo free energy of S are quite similar to those in 
Ref. [|13j - we will therefore be quite brief. We impose the length constraint on the n a field 
by a Lagrange multiplier; at N = oo this Lagrange multiplier is frozen at its saddle-point 
value and gives a 'mass' m to the n a field. The value of m is determined by solving the 
saddle point equation, which at t — t c is (using units in which % = kg = c = I) 



r d 2 k 

2—1 j A„2 



1 - 2p 



2p 



47r 2 y k 2 + u 2 + m 2 k 2 + (uj n — igH) 2 + m 2 



d 3 q 1 



(4.16) 



It is easy to check that m = is a solution at T = H = 0, confirming that the system is 
indeed at t = t c . As in Ref. [0] it can be shown that the leading term in the solution for m is 
independent of the nature of the ultra-violet cutoff. Evaluating the frequency summations, 
and a subsequent momentum integration we find that (^4.16 ) reduces to 



1 - 2p) log (l - e~ e ) + p log (l - e' e ~ r ) + p log (l - e 



-0+r 



where 



knT 



T 



e 

~2 



(4.17) 



(4.18) 



The Eqn. ( [4.17|) implicitly determines B as a function only of r. 
The N = oo result for the free energy is 



Jf = | E / ((1 - 2p) log (k 2 +u£ + m 2 ) + 2p\og (k 2 + {u n - igHf + m 2 )) - ^ 



29 
(4.19) 



We evaluate (|4.19|) using the methods of Ref. |TJ[ and find (after reinserting factors of k B , 
h, c) 



(4.20) 



with 



G 3 1 



n(0 = J^~^J e ydy[i l - 2P) log (1 - e^j + plog (l - e"^) + plog (l - e ^+ r 



(4.21) 



where B is also a function of r specified by ( |4.17| ). Exact evaluation of the integrals in Q(r) 
is not possible, but we have obtained the following asymptotic results 
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fi(r) 



2C(3) v^P, 
+ log 



57T 



2tt 



12tt 



— r + 
12 



'y/5 + r 
K(3)' 



2tt 



r -> 



r — > oo 



(4.22) 



where ( denotes the Reimann zeta function. In obtaining the above result we have used the 



non-trivial polylogarithm identities discussed in Ref. [29 



Results for the scaling functions for the specific heat and magnetization now follow as 
before and are plotted in Fig [7]. The Wilson ratio (Eqn ( |3.8|) ) can be obtained by taking 
the appropriate ratio, and the results are shown in Fig || The scaling function has the 
asymptotic limits 



f 5y/Ep / yg+ 1 ' 



n w (r) = < 



12C(3) 



2p7T 5 



r -> 



r — > oo 



(4.23) 



V. PHASE TRANSITIONS IN QUANTUM FERROMAGNETS 

We now consider the application of the ideas of this paper to one of the very first models of 
quantum phase transitions that was considered by Hertz |IJ : the zero temperature transition 
from ferromagnet to a paramagnet in an itinerant Fermi gas. The order parameter for this 
transition is clearly the local magnetization density, m a {x 1 r) (a = 1,2,3). This transition 
is special in that m a has a dual role - it is also the conserved charge density associated with 
global spin rotation invariance. 

The order parameter susceptibility 

x(x,t) = (m a (x,r) ■ m o (0,0)) (5.1) 

is expected to satisfy the following homogeneity relationship at the quantum fixed point 

x'(x',r')=s d+z - 2 ^ X (x,T) (5.2) 

This relationship defines the value of the critical exponent rj. The scaling dimension of m a 
is then immediately fixed at (d + z — 2 + rj)/2. However, m a is a conserved charge density, 
and below the upper critical dimension its scaling dimension must be precisely d. Equating 
the two scaling dimensions we get one of our main results 

z = d + 2-r] (5.3) 

Thus the three independent exponents z, rj, v have been reduced for the paramagnet- 
ferromagnet transition to just two - the values of z and rj are no longer independent. 

It is not difficult to see that Hertz's simple paramagnon model in fact violates the expo- 
nent equality ( |5.3j ) below its upper critical dimension. In his model 

z = 3 + 0(e 2 ) V = + O(e 2 ) (5.4) 
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where e = 1 — d is the deviation from the upper critical dimension. It is clear that these 
relationships are inconsistent with (|5.3f) at order e. 

We believe that this discrepancy can be traced to a more basic difficulty with Hertz's 
effective action: the incomplete treatment of the Ward identity ( |2.22|) associated with total 
spin conservation. In the vector-formulation (equivalent to our m a ) of the effective action 
there clearly must be cubic terms present if (|2.22|) is satisfied; such terms our absent in 
Hertz's treatment. Hertz also has a scalar-field formulation in which no such cubic terms 
will arise; however the full symmetry of the effective action is then hidden, and one cannot 
expect a proper treatment of the critical phenomena. 

A complete analysis of this problem clearly requires a more detailed consideration of 
the effective action of paramagnons, including the effect of cubic paramagnon vertex T. A 
preliminary analysis along these lines suggests that the identification of the upper critical 
dimension of d = 1 is incorrect. 



VI. CONCLUSIONS 

This paper has discussed the theory of the non-linear response to an external field, H, 
of a bulk quantum system in the finite temperature quantum- critical [TT,^fT2,r3] region of a 



zero temperature, second-order phase transition. In particular, we have considered the case 
where H couples to a conserved charge. For such a field, we obtained the general result that 
the scaling dimension of H is equal to that of ksT, even in the absence of hyperscaling. This 
result is encapsulated in the scaling form (|3.3|) . We also introduced a generalized Wilson 
ratio ( |3.8| ) associated with the non-linear response, and argued that it was a fully universal 
function of H/ k%I '. These principles were illustrated by calculations on some model systems. 

Tsvelik and collaborators have also recently studied the non-linear field dependence 
of the thermodynamics of heavy-fermion alloys. However they did not consider the special 
consequences of having a total conserved spin. The experimental data appear to indicate 
that the scaling dimension of H is unequal to that of k^T J7|. Using our results we may 
then conclude that any theory with a conserved total spin (some of the speculative proposals 
in Ref. || have a conserved spin) cannot explain the data. Spin-orbit scattering from the 
impurity sites must be included in an essential way in the final theory. 
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FIGURES 



FIG. 1. The universal scaling function Qm and Qc for the magnetization and specific heat (Eqn 
[)) for the Bhatt-Lee model with a = 0.6. The functions are obtained by taking appropriate 



derivatives of (4.8). The coordinate r = gH/(kBT). 



FIG. 2. The universal scaling function f° r the Wilson ratio (Eqn ( |3.8| )) for the Bhatt-Lee 
model with a = 0.6. 

FIG. 3. Ground states of the 2 + 1 dimensional O(N) sigma model in a magnetic field H 
described by the action ( 4.14| ). We have specialized to N = 3, p = 1/3. The coupling t measures 



the strength of the quantum fluctuations. The H = critical point at t = t c has the dynamic critical 
exponent z = 1. The line separating the quantum-disordered and XY ordered phases represents 
second-order transitions with z = 2. This phase boundary approaches H = as H ~ (t — t^) v 
where v is correlation length exponent of the classical Heisenberg ferromagnet in three dimensions. 

FIG. 4. Finite temperature properties of the model of Fig |3| for t < t c . There is a Koster- 
litz-Thouless transition at Tkt separating a phase with algebraic XY order from complete disor- 



der. The dependence fo Tkt on H at small H can be deduced from the results of Ref. [27]; here 
p s is the fully renormalized spin stiffness of the Heisenberg order at T = 0, H = 0. 

FIG. 5. Finite temperature properties of the model of Fig ^ for t = t c . The number K, is 
universal. The small T,H properties are described by the z = 1 critical point at T = 0, H = 0, 
t = t c and obey the scaling form (|3.3|). The scaling function Q(r) has a singularity at r = K,. 



FIG. 6. Finite temperature properties of the model of Fig || for t > t c . The physics of this 
phase diagram is discussed in some detail in Ref. [26|. The dashed line represents a crossover, 



while the full line is a Kosterlitz-Thouless transition. The functional form of Tkt is deduced from 
Ref. p|. 



FIG. 7. The universal scaling function and f^c for the magnetization and specific heat 
(Eqn (|3.7|) ) for the O(N) model in a field (|4.14| ) at t = t c . The results are obtained in the large N 
limit and plotted for N = 3, p = 1/3. The scaling functions are obtained by taking appropriate 
derivatives of (4.17 J4~2~I ). The coordinate r = gH/(kBT). The actual scaling function for pN = 1 
will have a weak singularity at r = 1C (corresponding to the Kosterlitz-Thouless transition of Fig ||) 
which does not appear in the large N calculation. 



FIG. 8. As in Fig [/], but with the results for the scaling function £lw for the fully universal 
Wilson ratio (Eqn ( |3.8|) ). The asymptotic limits of are given in Eqn ( p3|) . 
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